    p_rgh.storePrevIter();
    rho.storePrevIter();
